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Abstract 



Scalar field theory with an asymmetric potential is studied at zero temperature and 
high-temperature for theory with both ip and symmetry breaking. The equa- 
tions of motion are solved numerically to obtain 0(4) symmetric and 0(3) cylindrical 
symmetric bounce solutions. These solutions control the rates for tunneling from the 
false vacuum to the true vacuum by bubble formation. The range of validity of the 
thin-wall approximation (TWA) is investigated. An analytical solution for the bounce 
is presented, which reproduces the action in the thin-wall as well as the thick-wall 
q . limits. 
00 

I. INTRODUCTION 

The problem of decay of a metastable state via quantum tunneling has important applications 
in many branches of physics, from condensed matter to particle physics and cosmology. The 
tunneling is not a perturbative effect. In the semi-classical approximation, the decay rate 
^ ; per unit volume is given by an expression of the form 

H 

F»\ r = A e~ SB , (1) 

where Se is the Euclidean action for the bounce: the classical solution of the equation of 
motion with appropriate boundary conditions. The bounce has turning points at the config- 
urations at which the system enters and exits the potential barrier, and analytic continuation 
to Lorentzian time at the exit point gives us the configuration of the system at that point 
and its subsequent evolution. The solution of the equation of motion looks like a bubble in 
four dimensional Euclidean space with radius R and thickness proportional to the coefficient 
of the symmetry breaking term in the potential. When there are more than one solution 
satisfying the boundary conditions, the one with the lowest Se dominates Eq. (p. The 
prefactor A comes from Gaussian functional integration over small fluctuations around the 
bounce. The zero-temperature formalism is well-developed |l|, [|. In particular, it has 
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been proved rigorously that the least action is given by the bounce which is 0(4) invariant 

Linde [ffl] extended the formalism to finite temperatures. He suggested that at tempera- 
tures much smaller than the inverse radius of the bubble at zero-temperature, the bounces 
are periodic in the Euclidean time r direction and widely separated. Beyond this temper- 
ature they start merging into one another producing what is known as "wiggly cylinder" 
solutions. As one keeps increasing the temperature these wiggles smoothly straighten out, 
and the solution goes into an 0(3) invariant cylinder (independent of Euclidean time r) 
solution that dominates the thermal activation regime. 

The model field theory considered by most authors is p A theory with a symmetry-breaking 
term proportional to (p. However, it was shown || that the leading temperature-dependent 
correction to the tree-level scalar potential is proportional to p> 3 . Thus it becomes of interest 
to obtain bounce solutions to p> A theory with p 3 symmetry-breaking. An interesting special 
case is the so called thin- wall approximation (TWA), when the bubble radius R is much 
larger than the thickness of the bubble wall. In this limit, there is an analytical formula for 
Se in terms of the wall surface energy, and the details of the field theory are unimportant. 
However, it would be nice to also have an analytical interpolating form for the solution itself. 
In the TWA limit, we expect the bounce solutions in the two theories to be same, but for 
sufficiently large coupling constants for the asymmetric terms we expect differences. Also, 
it is not clear a priori what the limit of validity of the TWA is. 

Some of the above questions have been addressed by Adams . By suitably scaling and 
shifting the field <p as well as the potential U((p), he puts a general quatric potential in the 
form 

U(<p) = \p A - p? + (2) 

where < 5 < 2. He then computes the 0(4) and 0(3) bounce action for various values of 5. 
The thin- wall limit corresponds to 5 = 2, and the TWA result is recovered in that case. This 
throws up the fresh question: how does the difference between <p and p 3 symmetry-breaking 
get reflected in the values of 5 ? What is a thin wall in each of the two cases? 

In this paper we address the above issues. We obtain accurate numerical solutions for the 
zero-temperature and high-temperature bounces for p 4 theory with both p and ip 3 symmetry- 
breaking. We compute the actions in each case, and find that, for a modest value of the 
asymmetric coupling /(= 0.25), the action given by TWA formula agrees to within 9% 
with that obtained from the numerical solution for p breaking. We also relate these values 
to the 5 parameter of Adams However, the agreement is considerably better for p 3 
breaking than p> breaking. In the former theory, we verify that as / is reduced the error in 
the TWA formula goes to zero. We propose another criterion for the goodness of TWA, in 
terms of the temperature T* at which the actions of the 0(4) and 0(3) solutions become 
equal. A numerical investigation shows that the TWA holds up to / ~ 0.30. Finally, we 
present an analytical solution which satisfies the equation of motion with parameters fixed by 
demanding stationary action. This reproduces TWA results very well and, in the thick-wall 
limit, is in fairly good agreement with the numerical results. 

In Sec. II we review bubble formation theory at zero and finite temperature, and discuss 
the thin- wall approximation (TWA). In Sec. Ill we present our numerical results for 0(4) 
and 0(3) cases for p> and p> 3 type symmetry breaking as well as our investigation of the 
region of validity of TWA. In Sec. IV the analytical solution is presented in the limits of a 
thin wall as well as a thick wall. Finally, Sec. V contains some concluding remarks. 
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II. BUBBLE FORMATION 



Let us consider a scalar field theory with a Lagrangian density 

C{ip) = \{d^?-U{<f), (3) 

where the potential U(p) has two minima at p_ (false vacuum) and p> + (true vacuum). 

In the semi-classical approximation the barrier tunneling leads to the appearance of 
bubbles of a new phase with p = p + as classical solutions in Euclidean space (i.e., imaginary 
time r). To calculate the probability of such a process in quantum field theory at zero 
temperature, one should first solve the Euclidean equation of motion : 

with the boundary condition p — > v 9 - as x 2 + r 2 — > oo , where r is the imaginary time. The 
probability of tunneling per unit time per unit volume is given by 

r = A e" 5BM , (5) 

where 5^[<^] is the Euclidean action corresponding to the solution of Eq. (|j) and given by 
the following expression : 



SeVp] = J d 4 x 



(6) 



It is sufficient to restrict ourselves to the 0(4) symmetric solution p(x 2 + r 2 ) , since it is 
this solution that provides the minimum of the action S^kp] ||. In this case Eq. @ takes 
the simpler form 

fp + 3 dp = dU(p) 
dp 2 p dp dp ' 

where p = yjx 2 + r 2 , with boundary conditions 

dp 

p — > p- as p — > oo, — — = at p = 0. (8) 

dp 

We denote the action of this solution by S , 4 . 

Now let us consider the finite temperature case. Following ||], in order to extend the above 
mentioned results to high temperature, T ^ 0, it is sufficient to remember that quantum 
statistics of bosons (fermions) at T ^ is equivalent to quantum field theory in the Euclidean 
space-time, periodic (anti-periodic) in the "time" direction with period (3 = T -1 . One 
should use the T-dependent effective potential U((p, T) instead of the zero-temperature one 
U(<p) = U(p>, 0). Instead of looking for 0(4)-symmetric solution of Eq. ([|), one should look 
for 0(3)-symmetric (with respect to spatial coordinates) solutions, periodic in the "time" 
direction with period (3 = T~ x . At sufficiently large temperature compared to the inverse of 
the bubble radius R at T = 0, the solution is a cylinder whose spatial cross section is the 0(3)- 
symmetric bubble of new radius R(T). In this case, in the calculation of the action Se(p) the 
integration over r is reduced simply to multiplication by T -1 , i.e., Se[p] = T^ 1 S$[(p\. Here 
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SeI'p] is the four-dimensional action and Ss[ip] is the three-dimensional action corresponding 
to the 0(3)-symmetric bubble and given by : 



Sa[(p] 



d 3 r 



-(V^) 2 + f%,T) 



To calculate Sz(ip) it is necessary to solve the equation 

cPtp 2 dip dU(<p,T) 



dr 2 ^ r dr 



dip 



with boundary conditions 

<f — > <p_ as r 



oo, 



dtp 
dr 



at r = 0. 



(9) 



(10) 



where r = ^/S 2 . The complete expression for the probability of tunneling per unit time per 
unit volume in the high-temperature limit (T » R~ l ) is obtained in analogy to the one 
used in M and is given by: 

T(T) = A(T) e~ s ^ T]/T . (12) 

In the theory of bubble formation , the interesting quantity to calculate is the probability of 
decay between ip = ip_ and tp = tp + which are the two minima of U(ip). There is an interesting 
case (in the sense that the action can be calculated analytically) when U{<p> + ) — U{ipJ) = e 
is much smaller than the height of the barrier. This is known as the thin-wall approximation 
(TWA). At T = 0, in the TWA limit, the action S4 of the 0(4)-symmetric bubble is equal 
to 



Sa 



2tt< 



dp p 3 



= -^en 2 R 4 + 27r 2 R 3 S 1 . 
Here Si is the bubble wall surface energy (surface tension), given by 



Si 



dp 



dp 



(13) 



(14) 



and the integral should be calculated in the limit e —* . 

The bubble radius R is calculated by minimizing S4 with respect to R and this gives us 



R 



3Si 



whence it follows that 



5, 



27n 2 Sf 



The condition for the applicability of TWA is the following 



3Si 



» (d 2 f%-)/dy? 



2X-1/2 



(15) 
(16) 

(17) 



where (d 2 U (<£>-)/ dip 2 ) 1//2 is simply the order of magnitude of the bubble wall thickness. The 
results presented above were obtained by Coleman H. 
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These results can be easily extended to the case T » R 1 [||. To this end it is sufficient 
to take into account that 



An / dr r 2 
Jo 

4 



^ettR(T) 3 + 4 7 r J R(T) 2 5 1 (T), (18) 

where S\(T) is the bubble wall surface energy (surface tension) at finite temperature and is 
given by: 



S 1 (T) 



dr 



,d(p 2 
dr > 



+ U((p,T) 



(19) 



As before, the integral should be calculated in the limit e — > 0. 

The bubble radius R{T) is calculated by minimizing S3 with respect to R(T) and this 
gives us 

25i (T) 



R(T) 



whence it follows that 



lQnSfjT) 
3^2 ' 



The condition for the applicability of TWA is the following 



2Si(T) 



>> {d 2 U{ip^,T)/dip 



2N-1/2 



(20) 
(21) 

(22) 



where (d 2 U((p_,T)/d{p 2 ) x l 2 is the order of magnitude of the bubble wall thickness at high 
temperature. 



III. NUMERICAL RESULTS 

We start with the Euclidean action at T = 0, 

Se[p] = J d A x 
If we have 0(4) symmetry, Eq. ( p3|) reduces to 

S 4 = 2tt 



- - 1 dp p 3 



(23) 



(24) 



We compute the action for different values of the parameter in the symmetry-breaking term 
in the potential U((p). 

First we consider the following potential : 



U(<p) = ^tf - I* 2 ) 2 - F<p, 



(25) 



where A is the coupling constant and p is the mass of the scalar field. We perform the 
following rescaling in the action 



if -> if/p, p -> pyf\fj?, f 



Xp 3 



(26) 



to get 



271-2 

S4 = — / dp p 3 



2 ( ^ } + 2^ 



(27) 



The only adjustable parameter in the Lagrangian is /, so by covering the whole range we 
should be covering all relevant cases. 
The equation of motion is now 



d?tp 3 dip 



+ = 2{<p a - <p) - f, 



dp 2 p dp 

and the boundary conditions are the usual ones : 



(28) 



(29) 



dip 

ip — > ip- as p — > 00, — — = at p = 0. 

dp 

By solving Eq. (|2q) numerically for different values of /, substituting the solution in Eq. 
and integrating, we obtain the action for each value of /. 

At high temperature, we look for the 0(3) symmetric solution with cylindrical symmetry. 
Then Eq. (ffi) takes the form 



S* =4tt 



dr r 



l ( d<Ps2 
2 { dr ) 



+ U(<p) 



(30) 



Using the same potential (|25|) and the rescaling (|26|) , we get 



S, 



4np 



dr r 



1 .dip s 



V\ Jo 

The equation of motion is then 

d 2 ip 2 dip 
dr 2 r dr 



(31) 



2(^ 3 -¥»)-/, 



and the boundary conditions are the usual ones 

<p — > ip_ as r — > 00, 



dip 
dr 



at r = 0. 



(32) 



(33) 



Again, we solve Eq. fl32|) numerically for different values of /, substitute the solution in 
Eq. ( |3lD and integrate to obtain the action for each /. 

Table I shows our numerical results at T = and high temperature for different values 
of the asymmetry parameter /. 

Table I 



f 


S4 (Numerical) 


53 (Numerical) 


^(Analytical) 


S3 (Analytical) 


0.25 


3064.70 


143.80 


3368.80 


158.90 


0.55 


171.30 


18.70 






0.75 


13.20 


1.60 
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As we discussed in section II, for small values of / we can use the TWA formula for computing 
the action. From Eq. fll~4] ) 



Si 



dr 



o 



dr 



~ J dlf yf (if 2 

4 



The radius is given by 



R 



35i 



(34) 
(35) 



where e — 2/ (see @). For / = 0.25, we have R = 8 and the value of the action is (see 
Eq. ©) 

S 4 = 3368.80. (36) 

Comparing this analytic value with the numerical value for / = 0.25, we get an error equal 
9%. 

At high temperature, again Si = 4/3. The value of R(T) = 16/3 and the action is (see 
Eq. (P) 

S 3 (T) = 158.90. (37) 

Comparing this analytic value with the numerical value for / = 0.25, we get an error equal 
9%. Thus even for / as small as 0.25 the TWA formula for the action does not give very 
accurate results. Obviously, there is no point in comparing numerical results obtained for 
higher values of / with the TWA formula. 

Now we repeat the above calculations for the following potential 



c%) = ^ 2 - p 2 f - iV. 



(38) 



We scale ip> and p as before, with the dimensionless coupling / defined by / = Fj (Xp), solve 
the equation of motion and calculate the action. Table II shows our numerical results for 
different values of /. 

Table II 



f 


S4 (Numerical) 


S3 (Numerical) 


S 4 (Analytical) 


S 3 (Analytical) 


0.25 


3446.10 


161.80 


3368.80 


158.90 


0.55 


349.50 


35.40 






0.75 


139.20 


19.90 







We now check the TWA formula for / = 0.25. From table II we see that the error in S4 
is only 2.3%, while the error in S3 is 1%. Thus we see that for / = 0.25, the TWA holds 
much better for symmetry breaking with a (p 3 term than with a (p term. 

Notice that the TWA formula is independent of the symmetry-breaking term in the 
potential. Moreover, the action for large values of / for the (p 3 symmetry-breaking term is 
larger than that for the (p one for both the 0(4) and the 0(3) cases. 
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To test our numerical method (we have used Hamming's modified predictor-corrector 
method for solving the equation of motion), we have calculated the action for small values 
of the symmetry breaking parameter / in the potential (|38|) and compared it with the TWA 
formula. In Fig. (1), we plot the percentage error in the TWA formula as a function of /. 
The dots represent our results while the solid line shows a fit to the data. We see that the 
error decreases for small /, as expected, and approaches zero as / — > 0. 

Table III shows the relationship between the values of 5 in Eq. (0) and / for ip and (p 3 
perturbation. 



Table III 


/ 


6 


ip case ip 6 case 


0.0 


2.0 2.0 


0.25 


1.85 1.85 


0.55 


1.50 1.62 


0.75 


0.64 1.46 


0.77 


0.0 1.44 


0.80 


1.41 


1.0 


1.25 


5.0 


0.15 


10.0 


0.03 



As we can see from the table, as / goes from to 0.77, 5 goes from 2 to for the (p case 
(for / > 0.77 the maximum disappears) and from 2 to 1.44 for the tp 3 case. Also we can see 
that 5 goes to for the p 3 case only at large values of / (here the maximum disappears only 
asymptotically). 

As already mentioned, at zero temperature the 0(4) symmetric solution has the lowest 
value of Se, i.e., Sg = At high temperature, we have S& = S 3 /T. At intermediate 
temperatures other solutions exist. In the TWA, however, it has been shown || that all other 
solutions have higher Euclidean action. This corresponds to a first order phase transition 
from quantum tunneling at low temperature to thermal hopping at high temperatures. The 
transition temperature T* is given by equating S4 with S3/T, i.e., 



S3 
04 



If the surface tension Si is temperature independent, we have 

27n 2 Sf 



S 4 -- 
S 3 



2s 3 

16ttS 3 
~3^~ 



Dividing Eq. (g0|) by Eq. (gl]) and putting e = 2f (see 0) we get 

T„ = C*f 

where 



C 



64 

8171-S! 



(39) 

(40) 
(41) 

(42) 
(43) 
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Thus we see that, in the TWA, T* increases linearly with /. We test this by computing 
S3/S4 from our numerical solutions at different values of /. Fig. 2 shows our results for 
the potential of Eq. (|38|). We see that, for / < 0.3, there is very good agreement with the 
predicted linear dependence. This also confirms that, in the domain of validity of the TWA, 
the surface tension S±(T) is independent of T. Beyond / ~ 0.3 in our dimensionless units, 
there is a systematic deviation from linearity. Thus we can say that, for values of / larger 
than this, the wall thickness becomes important. 



IV. ANALYTIC SOLUTION IN TWA LIMIT 

We use the following potential (see ||) to calculate the action analytically in two extreme 
limits: the thin-wall and thick-wall limits. 

U{<p) = \p' - ^ 3 + \^ 2 . (44) 

Thin-wall limit : 5 — > 2 

We find that an analytic solution for the bounce of the form of a Fermi function: 

1 

V = e ( P 2 -R 2 )/A 2 + p ( 45 ) 

where p = yjx 2 + r 2 , R is the radius of the bubble and A its width, acts like a bounce in the 
TWA and leads to the correct value for the action S4. 

Here the false minimum of the potential is at <p = and the true minimum lies between 
2 (for 5 = 2) and 3 (for 5 = 0). The parameter 7 is approximately equal to true minimum 
in the TWA. The bounce has values ip = 7 at p = and at p — > 00. The boundary 
conditions (§) are satisfied by Eq. (|45|) . 

We discuss the zero-temperature action S4. To evaluate 7, R, and A, we substitute the 
ansatz (gD in Eq. (0) : 



d ip 3 dip 



Then the left-hand side (L.H.S.) and the right-hand side (R.H.S.) are respectively 
LHS 78p 2 /A 4 + 7 (-12p 2 /A 4 + 8/A 2 ) 



( e (p2-fl2)/A2 + ^ 3 ' (g(p2-/?2)/A2 + ^ 



2 



7 (4p 2 /A 4 - 8/A 2 ) 

+ e (p 2 -R 2 )/A 2 + 1 ' { l) 



7<5 37 2 

R H S - = e{p 2_ H 2 )/A 2 ~ j - ( e (p2_R2 )/A 2 + ^ 2 

7 3 

+ ( e (p 2 -i? 2 )/A 2 + fj3" ( 48 ) 

In the TWA, the solution is constant except in a narrow region near the wall at p = R. So, 
we replace in Eq. (|47|) 

8R 2 1 
8p 2 /A 4 by -^(l- a A 2 /R 2 ) in the (e(p2 _ R2)/A2 - 1)3 term, (49) 
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VA 2 -12p 2 /A 4 by 



12R 2 



1 - bA 2 /R 2 ) in the 



4R 



( e (p 2 -R 2 )/A2 + iy 



4p 2 /A 4 -8/A 2 by ^(l-dA 2 /R 2 ) in the 



1 



g(p2-R2)/A2 _|_ X 



term, (50) 
term, (51) 



where a, 6 and d are parameters to be determined later. 

Comparing Eq. (g7p with Eq. (g§) in the range R 2 (l-A 2 /R 2 ) = R 2 -A 2 < p 2 < R 2 +A 2 
R 2 (l + A 2 /R 2 ) where p 2 ~ R 2 as A 2 /i? 2 << 1 , we have : 



7 2 


R 2 


8 


~ A 4< 


7 


R 2 


4 


~ A 4 


6 


i? 2 


4 


~ A 4 



, 2 / D 2> 



i 2 / ^,2^ 



(52) 



(i - rfA 2 /i? 2 ; 

We can now evaluate the zero-temperature action S4 : 



Sa = 2tt 2 



dp p 2 



(53) 



Substituting Eq. ( f4"5| ) in Eq. ( p^) and integrating we get 
S, 



2w 2 7 2 R 4 



6A 2 



,T 2 sA\ 5, a2/ o 7T 2 A 4 

l + (--2)— ) + -(!- 2A 2 /it> 2 + 



3 i? 4 ' 



+2(1 _ 3A 2 //? 2 + (y + l)A 4 /i? 4 ; 



16 V 



11A 2 ~ 2 



,7T 



+ (— + 2)A i /R 4 ' 

3R 2 V 3 ; 7 ' 



(54) 



We now determine the parameters a, 6, and d by demanding dS^/dR 2 = dS^/dA 2 = dS^/dj = 
0. Differentiating Eq. fl54 ) and using Eq. (|52|) , we find that to leading order in A 2 /R 2 , 



4b - 2a - 2d + 1 = 0, 

36 - 2a - d = 0, (55) 
36 - lla/6 - d = 0, 

which leads to a = 0, 6 = 1/2 and d = 3/2. Using Eq. (|52j), we can rewrite Eq. fl54]) as : 



5 4 = 2tt 



,8i? 6 



A 2 



(1/3 - d/2 - a/2 + 6) + — (d-36+ lla/6) 

R 2 



(56) 



where the coefficient A^ evaluated by the usual methods of statistical mechanics for the 
Fermi function vanishes. This gives 
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i? 6 



(57) 



The quantities 7, R and A are determined from Eq. ([52]) using the values of a, 6, and d. So 
we have 

T 2 -27^+2*_ = 0, (58) 

which gives 

A' = ^ = (59) 
7 Z — 27 7 — 

with 7 given by Eq. (|58|) . We have then, for 5 = 1.9, 7 = 2.1, which implies that R 2 /A 2 = 

-^r + b = 11. Thus we have 

7—0 

47T 2 

^4 = — (ll) 3 , (60) 

while the action from the TWA formula is ( see 0) S TW = ^f(10) 3 for 6 = 1.9. The 
departure from TWA, S4/ Stw — 1-33, is in agreement with Ref. |J. The expressions seem 
certainly valid for values of S in the range 2.0 to 1.8 . 



Thick-wall limit: 5 — > 

The form of the bounce in Eq. ( [45]) suggests that the thick wall limit, which would 
correspond to small values of R 2 /A 2 , would be obtained by approximating the Fermi function 
by the Maxwell-Boltzmann function, which leads to a Gaussian: 



(p — 7e 

The action for this form of bounce is found to be 



P 2 /A 2 



C 2 2 a 4 

04 = 7T 7 A 



1 5 7 7 2 
+ ;r - VT + VTT 



2A 2 8 9 64 



(61) 



(62) 



Eqs. ([52] ) then reduce to 



r 



a 

A 2 



2 

4 



A 2 



4 



d 

A 2 



(63) 



Note that in this case 7 1, so 7 2 is negligible. 

The values of b and d are again obtained by demanding dS^/dA 2 = dS^/dj = 0. This 
gives b = —9/8, d = —1/2, giving 



A 



9 A 



This yields the action 



4vr 2 R 2 
S, = — (1-9)5 + 0(^). 



The ratio of the action to the TWA value is 

S 4 



R 1 



TW 



1.95(2 -Sf 



(64) 



(65) 



(66) 



For 5 = 0.1, i?4 = 1.31, which agrees with Adams' result. 
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Thus, the form of the bounce given by Eq. fl5"2]) seems valid over the whole range of 5 ( 
from to 2 ), and in the two extreme limits is amenable to analytic calculations. 

The numerical plots of the function ip in the two extreme limits are shown in Figures 3 
and 4. In Fig. 3 we compare our numerical result with the analytic one for 5 = 1.8. From the 
figure we see that the Fermi function agrees very well with our numerical results. Similarly, 
Fig. 4 shows the Gaussian function compared with our numerical results for 5 = 0.2. We 
find that 7 = 0.55 numerically and 0.45 analytically. This discrepancy in the values of 7 is 
expected, because the terms neglected in Eq. ( |55]) are of order (pr). 

The analysis given above can be carried out for the 0(3) symmetric bounce and we 
find that this gives equally good results. We anticipate that this approach will be of great 
use is discussing the extension to arbitrary temperature. This will be discussed in a future 
publication. 

V. CONCLUSIONS 

We have obtained accurate numerical solutions for the zero-temperature and high-temperature 
bounces for both (p and ip 3 symmetry-breaking. We compute the actions in each case and 
find that, for a modest value of the asymmetric coupling / = 0.25, the action given by the 
TWA formula agrees to within 9% with that obtained from the numerical solution for (p 
breaking while for ip 3 the agreement is within 2.3%. Hence, the agreement is considerably 
better for (p 3 breaking than ip breaking. At high temperatures, the conclusion is qualitatively 
similar. 

We have checked our numerical method by comparing the action obtained numerically 
with the one obtained from the TWA formula. Very good agreement is obtained as we go to 
small values of /. We also verify that as / is reduced the error in the TWA formula goes to 
zero. We propose an alternative criterion for the goodness of TWA, in terms of the relation 
between / and the temperature T* at which the actions of the 0(4) and 0(3) solutions 
become equal. A numerical investigation shows that TWA holds up to / ~ 0.30. Finally, 
we present an analytical solution which satisfies the equation of motion in an approximate 
sense in two limiting cases. The first of these reproduces the leading corrections to the TWA 
results very well. The second is applicable for the opposite case of a very thick wall. This 
gives us insights into the nature of the bounce solutions for various values of S going from 
thin to thick walls. 

Our work overlaps to some extent with that of Adams ||. Like us, he had computed the 
bounce action for both 0(4) and 0(3) symmetry. Using the parametrization of Eq. (|^), he 
was able to map the entire range of wall thickness into the interval < 5 < 2. However, 
while the form (0) is computationally very convenient, it obscures the fact that the free 
parameter actually measures the departure from pure y? 4 theory, and does not directly allow 
a comparison between tp and ip 3 symmetry-breaking. Our parametrization allows such a 
comparison. Of course, both theories can be mapped onto (|2|), but the S — f relationship is 
different, as illustrated in Table III. Moreover, Adams does not attempt to obtain approx- 
imate analytical solutions, which is done here for the first time. Also, he does not lay any 
emphasis on the — f relationship and indeed does not compute T*. 

Much of the work on inflationary models relies on the zero-temperature potential, so 
our results could be relevant for inflation [[|. They may also have some bearing on the 
formation of topological defects in a first order phase transition where authors consider 
zero-temperature potentials, see for example [§]. 
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So far, we have discussed the action only at zero temperature and high temperature. 
To obtain the bounce solution at intermediate temperatures, we have to solve a partial 
differential equation with periodic boundary conditions in the r direction. We have done 
so using a multigrid method as well as by extending the analytic bounce calculations. This 
work will be presented in a future publication. 
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Figure Caption 

FIG. 1. Error in the TWA formula as a function of /. The squares represent our results 
while the solid line shows a fit to the data. 

FIG. 2. Deviation of T* from the TWA limit. The straight line represents the TWA limit 
while the squares are our numerical results. 

FIG. 3. ip as a function of p. The solid line is the Fermi function while the dashed line 
is the numerical result. 

FIG. 4. ip as a function of p. The solid line is the Gaussian function while the dashed 
line is the numerical result. 
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